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ABSTRACT 


A brief account is given of a proposal by Dr. H. Schiff 
for a six-vector unitary field theory. The static approximation 
to the equations is found and is the basisfor the remaining part 
of the thesis. 

The nature of the spherically symmetric solutions of 
the static equation is discussed; most of the results being ta- 
ken from a paper by Finkelstein, et al., and it is shown that in 
addition to a continuum of solutions asymptotic to a there 
exist discrete solutions, asymptotic to Zero, which could repre- 
sent neutral bosons. 

The interaction potential energy between the discrete 
solutions as given by an approximate expression is discussed. 

The analytic interaction energy curves computed by Teshima are 
studied in detail for the existence of quantum mechanical and 
classical bound states. 

The static equation is shown to be the Euler equation 
of a certain Lagrangian. Several integral relationships satis- 
fied by the solutions are proved and on this basis a variational 
method which always gives upper bounds to the energy is formu- 
lated. 

The variational method is used to calculate spherically 
symmetric solutions and it is shown that these solutions do ap» 
proximate the spherically symmetric solutions obtained by numeri- 
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The research reported in this thesis is motivated by a 
propeial @f Drj4H. Sehift@r fer a unitary: sixz-vector field theory: 
By unitary field theory we mean that matter and its properties 
as well as the electromagnetic field are different states of the 
same physical entity, a "universal field". In Dr. Schiff's view 
the universal field is taken to be the electromagnetic field it- 
self, so that all matter is a manifestation of electromagnetic 
energy. A review of the previous attempts at such a theory, as 
well as some preliminary work on Dr. Schiff'ts theory will be 
found in a thesis by Teshima (1960). Since a number of Teshima'ts 
results will be used here, that work will be designated by T. 

As.a justification for a unitary field theory we quote 
the well known fact that quantum field theory is unable to pre- 
dictrthe mass or coupling constants for the particles introduced. 
Neither is there any limitation to the number of kinds of parti- 
cles. The hope is that these difficulties can be overcome by 
considering all particles as eigenstates of a universal field 
which is described by a universal field equation. Since only the 
universal field is present, all particles must arise from the 
interaction ef this field with itself. Therefore, the basic equa- 
tion must be non-linear. The universal field must also be a spin- 
or, otherwise we could not get solutions for the fermions, i.e., 
from spinors (representing fermions), one can always construct 


tensors (representing bosons), but the converse is not true. In- 
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vestigations along the above lines of reasoning have been carried 
out by Heisenburg (1957). Heisenburg's universal field has no 
classical counterpart. Let us point out again that in this attempt 
at a unified theory of matter the electromagnetic field is consi- 
dered as the universal entity. Since the electromagnetic field 

is a tensor, it is not reasonable to expect this theory to re- 
present fermions. . Furthermore, Dr. Schiff'ts theory at the present 
is classical; Heisenburg's is not. 

Several other theoreticians, notably Chew (1961), are 
convinced that the Pi-Pi interaction is fundamental for all ques- 
tions of strong interactions. Nevertheless, the dispersion rela- 
tions do not account for the dynamics of the interaction while, 
in principle, a non-linear theery would be able to account for 
the dynamics. 

We shall discuss the basic ideas of the theory only to 
a minimum extent, since we are primarily concerned with investiga- 
ting the solutions of the static field equation in this thesis. 


The fundamental equation is: 


= 
re: 2 — a : 
g 3x, - BA) + | ante = 0; Fi = ¢ =! (1 1) 
where Au is a four vector potential and a is defined by 
JA, OA | 
= a 1-2 
Pag ez, are (1-2) 


as in the usual electromagnetic theory, and g is a coupling con- 
stant with the dimensions of charge. The indices run from 1 to 4, 
with Xi» Xo» X3 = Me Vs 2s Xy = ict; the four potential is related 


to the ordinary three vector potential, A, and the scalar potential, 
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is the Klein-Gorden operator with a self interaction. 

Equation (1-1) is clearly non-linear in the four poten- 
tial. Only integer spin fields can be formed from tensors and 
Foo is a tensor with six independent components, so equation (1-1) 
is taken to represent bosons with integer spins including zero. 


The real and imaginary parts of equation (1-1) give 


the two sets of equations: 


re 2 
|v a -m - 2° (Ae ss °)| Foo = 9 (1-3) 
0 OA, 
2A. aes + oe Huy = 0 ae 


by Biv and not summing over uv. This will give: 


Oo “UD 


3 a 1-5 
| 4 F }-° (no sum; wp, v) (1-5) 


A spin zero field is described by a scalar function so as an ap- 
proximation to the spin zero case consider A= 0. By using equa- 
tion (1-2), the field equations will contain only the electro- 
static potential, o. if io O, then H = 0, hence the particles 


can have no magnetic moment. [From (1-5), with A= 0 


32 | (v0)? | = 0, 


and therefore the electrostatic potential is time independent. 


From equation (1-3) 
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Putting E = -Vo, 
2 
V 92 - m° + oe |e = 0, 


or, restoring hand ec, 


Ce . 
e_ mc” g 1 oe Ree 
L - 7 22 3? | 9 = Y= constant . (1-6) 


To get this equation in a different way, one might sim- 


ply postulate a spin zero field equation of the form of a Klein- 


Gorden equation with a non-linear self interaction: 


2 1 d<o ne em iS. 
ae ee a ee oo ee 
Get aay 1 e 


and consider the static limit. Let us point out that the sign 


of the non-linear interaction is important, since a positive 
sign as above will lead to particle like solutions, while a ne- 
gative sign will not. See N.V. Mitskevich (1956). 


To get a dimensionless equation make the substitutions: 


2 
aa _ me . = 
Dire ta ge ae ie SG ae 
1-7 
2 
me BG 
Ga and Py are the dimensionless potential and radial distance, 
respectively. Equation (1-6) becomes: 
7 © 
2 tae _g . 
b - 1 7s Pa Pq ~ ot ) (1-8) 


The subscript "dad" will be omitted in most of the following work. 
We shall only consider the special case of equation (1-8) with 
y= 0. The basic equation to be studied will then be: 

Vo -9+ee=0 . (1-9) 


If mo is to represent physical reality, » must be continuous and 
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bounded everywhere. 
In this work we shall investigate only those solutions 


with the following boundary conditions: 
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and for spherically symmetric solutions, (1-10) 


Any solution of equation (1-9) satisfying these boundary condi- 
Ti0Ons is called an eigensolution. 
We shall now proceed to discuss the nature of these ei- 


gensolutions. 
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II. Phase Plane Analysis 
of the Radial Equation 


Equation (1-9) was discussed in some detail by Finkel- 
stein et al. (1951). We shall reproduce their analysis below. 
We wish to investigate the solutions with spherical 


symmetry of the equation: 


Vp = - 5 o>. (1-9) 
We have: 
=, 55 (vr BB) = 9 -4e°. (2-1) 


as 


This equation has three trivial solutions: 

Go = Gonstany =f 
with 

Cai. Bae NS 4 (2-2) 
Let us find the behaviour of the solutions near these special 
points. To do this. sep 

gp =K+ U(r) (2-3) 
and regard U(r) small with respect to K. We can neglect the non- 


linear term obtained by substituting (2-3) into (2-1). We get: 


| S200 ,y Se oa = 
“3 ant ee ame TH) near K = 0, (2-4) 
leet PS CU owed (or pet m 
= So(r© 52) ~ = NB + U(r) - FI 33 + 9U(r)) (2-4B) 
= - P(r} near Ky eo V3 


In general the solution of 
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may be found by letting U = = then —= = - kf 
dr 
hence 
f = Ur = A sin(vk + q) fore >0 
eae for k< 0 


Then we see that: 
Near K, = T af 3 Where; kia aeS 
a (2-5A) 
the solution is trigonometric, 
Near K = 0 where k = = ] 
(2-5B) 
the solution is exponential. 

To investigate the solutions further let us now make a 
plot of what is happening in the phase plane. Let 9 and 9! = ge 
correspond to the "position" and "velocity" of a point in phase 
space. We are now thinking of the position r as "time" to obtain 
an analogy with the dynamical case. Then equation (2-1) corres- 


ponds to a non-conservative one dimensional motion, since r (the 


time) appears explicity. Rewriting equation (2-1): 


1 _ _ 2 ge 
S be O - Fe -te (2-6) 


oa: 
2 


(2-7) 


2 eee ee ee ee 2 (ae 
a Eo tie? nee? > (ar? 
When the motion is conservative the particle (of mass = 1) satis- 
fies. equations (2-6 and (2-7) with the right hand side set equal to 


zero. From equation (2-7) we see that the energy, E, for the con- 


servative motion is: 
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E = = (op! 
where V(r) = = 9 -s 2 
V(r) is the potential energy. The equilibrium points of the mo- 
tion are given by S = O and correspond to the special solutions 
po = O, aa: A point representative of the conservative motion 
moves along the lines of constant E in the phase plane. For given 
EK we can now make a plot of the conservative motion in the phase 
plane. This is shown in Figure 1. The curve for E = 0 is a figure 
eight through the origin; the curves E > O enclose both Ky and 
K-= Q, but the curves E < O enclose only one of the solutions KY 
oy f. 


The non-conservative motion, which corresponds to our 


actual problem may now be described. From equation (2-7) we see 


that $= = - E(S2)*. This term represents the dissipation of ener- 
sy with time. Therefore, since 24 is never positive the repre- 


dar 


sentative point of the actual motion must always move inward across 
the lines of constant E. Such a line must always end at either K, 
or the origin, no matter where it starts. In particular, if any 
solution curve gets into the shaded region, it must terminate on 
either KY or kK. 

The solution curve labled eo is designated as an "eigen- 
solution". These solutions are finite at the origin and approach 
zero as r—o. This follows from equation (2-5B). An eigensolu- 
tion is located by starting on the » axis and continuously in- 


creasing the initial value, 9(0). At first when 0 < 9(0) < V6 
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Figure 1 


A sketch of the phase plane for depicting solutions 
of radial symmetry. The curves beginning at kK and 
k_ are solutions asymptotic to V3 while the curve 


eo is an "eigensolution". 
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all solutions are certain to terminate on K,. The situation is 
still the same for (0) = k, which is larger than the critical 
value, V6. But for a value 9(0) = k_ the solution will terminate 
on K_. For some value 9(0): = e between k, and k_ the solution 
will terminate at the origin. By narrowing the interval Ky kK 
the eigensolution may be determined with arbitrary accuracy. 

The above analysis yields the following facts concerning 
the solutions of the equation V9 = -< > for radial symmetry. 


1. Every solution must behave at infinity either as 
Z) i snc as + a) ie 


=? 


or eS 


2, No solutions exist for which o'(0)=0 and (0) = 0 except 
the trivial one -o = 0. 

No eigensolution's solutions exist for '(0) and (0) < NGt 

4, If we require that o' be zero at the origin, then the solu- 
tions of the type a) form a continuous set, while those of 
type b) form a discrete set. 

5. The solutions forming the discrete set have zero, one, two, 
ete., nodes corresponding to the solution being the first, 
second, third, etc., eigensolution. 

The first three solutions of the type b) were computed 
by numerical integration of equation (2-1); the detailed discussion 
will be found in T. The solutions with no node, one node, and 
two nodes are designated by Po? Py and Po respectively, and the 


first two are reproduced here in Figures 2 and 3, 
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Figure 2 


The first elgensolution, Po» of spherical symmetry 
obtained by numerical integration of equation (1-9). 


Reproduced from Teshima, (1960). 
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Figure 3 


The second eigensolution, P42 of spherical symmetry 
obtained by numerical integration of equation (1-9). 


Reproduced from Teshima, (1960). 
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III. Bound States 


If the discrete solutions of equation (1-1) represent 
particles, it is important to compute the interaction between so- 
lutions, i.e., particles of different kinds. Note that this inter- 
action would be purely electromagnetic and in the static approxi- 
mation could be represented by an interaction potential energy 
curve. The difficulty in finding a proper expression for the 
energy of the field has been stressed in T. Since a Lagrangian 
for the theory has not been found, an "ad hoc" expression for the 
energy and the interaction energy of the field must be used. 


We assume for the energy the expression 


E = | poay Be Fa 
where p is the charge density, and satisfies 
-inp = Veo . (3-2) 


Upon converting the p and 9 to dimensionless quantities, using 
equations (1-7) 
2 
mc 
i = a 2 a | eavaae 
me“ 
is obtained; all of our energies will be in units of —- 


Extending (3-1) to two particles, we define the total 


energy as 


En = \(pz + Po) (%, + Po)dv . (3-3) 


When the systems are far apart we know that 
En = | p19,9¥ + \ pyepav , (3-4) 


Hence the interaction energy, Ey is given by the difference be- 
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tween (3-3) and (3-4), namely 
Es = [eyo av + Jeno, dv 


= 2 { pyepav = 2 | pp, av 
(It was shown in T., that [p,95av = { p99, av.) For an arbitrary se- 
paration, d, of the systems one and two, the interaction energy, 


Written in full isa: 


By(8) = 2 fo, (#ep(d - Fram 
or (3-5) 
o 2 foil? 7 d)o5(P) ar, ete 


(Note that Teshima's expression for the interaction energy differs 
from ours by the factor of two.) 

The possibility of forming bound states between parti- 
cles of the type Po and ?1 has been discussed in T. The interac- 
tion energy, Voy? for an arbitrary separation of these particles 
is plotted in Figure 4, This energy curve was obtained by fitting 
analytic functions to the numerical solutions of Po and ?1 and 
calculating the integral (3-5) An attempt was made in our work to 
compute the interaction by direct numerical integration using the 
numerical solutions Po and ?,- We had hoped to gain additional 
accuracy over Teshima'ts potential by eliminating the need to "fit" 
the numerical solutions with functions. The attempt failed, evi- 
dently because the numerical solutions were not computed for 
enough points, this caused the potential to undergo wild variations. 

If the bound state is formed classically, it is readily 


2 
seen from Figure 4 that the binding energy, B, will be 12.84 — 


{lemsa lie) bas ive 
vb -ogal + vbgo;a| = = 8 | 
wma 8 7 vige a |S: oii a 


-98 (Weatidas as 20% (. vga |= vii 4) tent ..T <a 
Naetens aottosietat ot .owt bae.sao emedeye oat to ee 


nek List on + oa 
| b("y - Dgohals = (B) 78. (ty. eae 17 a 
(ee) | ae 


ee ee aerate 
erstiih ygrecs noltorrednt arid aot aoteesigne atsmtdee? Bis: otoK) 
(owed ‘Yo aotos? ont Yd ‘dquo mort 
-liusq neswied sotste bnuod gitmrot Io YIEEtdteeog on? , | 
-ostedit ea? .? nf heseunath asad eatt 1? Bas op says odd To solo 
eslotisisq saodt ‘to nokisusqes’ Lala nal 107 “gov « 


gnisdt2 yd beatstdo sew svivo Yartens aid? ,# oxugtt rat pi a g 
bre .9 bns op To snoltutoz feo trons ont (98 enotvonuft ottytens eS 


we 


o¢ otow ao ot ebsm asw tqmetts nA 2-8) sone ‘entt Mie tal 
ed gatay neldstgednal sofas nae va nottostedat cd stuqnoo- 
fsnottibba ates ot beqod bad oW fp? bas 6@ anotsulos - 
"319" o¢ besa Ad satdsmtatis ya setsusiey: Camas? 20 
-lve ,beflst tqmetts enf - venottonut cttw enoktwtes | 
Jol bedugmoo dom stew enotsuioe homage 
 6foitsizrsy bitw ogrebriy od Latsaetog ont hanaee 2 
yltbssr at #2 wileoteeato’ hae ae osade & _ 


s Bi: 
om 
—- H8NSE od Eriw a rsa sata oe oi - al 


" 
eo Pe, 


Figure 4 


The black curve shows the interaction potential 


energy, Voy? 


second kinds. Reproduced from Teshima, (1960). 


between particles of the first and 


The square well approximation is shown by the 


dotted curve. 
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with a separation, r, of 1.QA. 
To solve the quantum mechanical problem for the binding 


energy a tractable form for the potential, Vo is needed. The 


-* 
Morse potential 


O 


V=V ai. as ~ ze" (-P0)/%0 (3-6) 


has been plotted on Figure 5 together with V The figure shows 


Om 
that this function is a good approximation to Vo. at least in 
the region of the well, where it is most important to have a good 
fit. Any further improvement in trying to fit an analytical curve 
tO Voi is unwarranted, since the existing Vo1 has been obtained 
by an approximation. 

The energy, Bo of the bound states for the potential 


(3-6) is found in Morse and Feshbach, Methods of Theoretical Physics 


(1953). © a; : 
i 
Pao Voth ae tp - Bre 
O ema, 
ho/a Ag 
for I= 0 and e ° Sa la: Wise any inveger ae NemV, - 3 


Converting these to the dimensionless form by using the relations 


(1-7): . - 5 
-_ Mie < 2 fe) I a 1 
B = me“ | —— + = |) se (2 + Se) oa (n + Se) 
n os eae 2" a2 2 
and (3-7) 
Ne i 
n< <8,\/oG ~ 2 


Following Teshima the energy of the 1st and end particles is ta- 


ken to be: (pages 53, 54 of T.) 


2 
mc ~ = 
Ey = =z (13.53); E, = 6.293 Eo (3-8) 
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pigure. 5 


The solid curve shows the interaction potential 


energy, V while the dashed curves are the 


Ol? 
Morse potential approximations. The parameters 
for the Morse potential curves are given in Ap- 


pendix D. 
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The number 13.53 and the ratio E/E, is not arbitrary since they 
are a direct result of the non-linearity of the field equations 
together with the fact that we are getting something similar to 
eigensolutions of a linear equation for the solutions Po and Pi: 
If..tmne,energy of the basic or first particle is to be. equal to 
me“, we must take the coupling constant to be 
G 2352. (3-9) 

Inserting numerical values, we find that for the three 
Morse potential curves plotted in Figure 5, n is any integer less 
than about 1/2. Therefore, n = 0, and there is only one bound 
state. 


Using (3-7), (3-8) and (3-9) the total energy, E of 


i? 
the bound system of particles Po and P41 can be calculated using: 


En = Ey - EL - BL (3-10) 


A table of values of Bo and Em is shown below for the 


three Morse potential curves plotted on Figure 5. 
Tapie 1 


Quantum Mechanical bound states between 
Po and P41 


| 
a 0.119 | 0.231 | 0.250 


a 


Energies are in units of me” with G = 13.53 


This does not exhaust the possibility of forming 
bound states between the solutions Po and Pi: As can easily be 


seen, if » is a solution to equation (1-9), then -p is also a 


{ gi 
7 es [ | 


ire ¢ ear 
- 7 


j < 7 2 pel me 
yout sonte ouabaniiads tor. ek Ng obtex ont bas | 


a 


ihe 


enotssups HLoLt ait ‘to ysimseatbesion ens ‘to pay | 
ot tslimts suntdvemos gnivdeg ste sw ted? tost sit wet ito: 
“(9 boas 5? anor tofoa srt 10 sotteups tsenti s to anol ul 
ot [sups ed ot at oLwkdtils tart? ro otesd sft to eat’ 
od ot tastenos antiques sit exad tasint oe 

(@-€) . €@.€£ = 2 


sexdd oft tol Jsdd, Daf ow ,eeutsv [sotvemun anidreanl 


7 


_ 


-) a 
se 


seal tozgetdt yas Bt o .2 esugli Gk beszolq sevawo setsneteg ¢ 
brived ono ¥YEhe et aisit bop .O = a .steTeted?® .S\f suods 1 


é 


~ 0% 5 
bea 
io «pi euptsns Isto sas (O-€)-fas (8-e) (T=€) snteU : _ 
‘antes bedsivoles ad dso 1? pris Q? seloliteg bod atic bayod | ais 

- oe 
ct 


i 
7 
> =) 


- (OL<e) i - 2+ =e 


ait 63 wolsd awoe al pe bas aa to aswisv to elds A a 


— 


.c siuala mo ottiite BavLuso Isttaetoq— oertoM 
[ sidsT 


agewied asiste baved L[solnsdooM mudaeup - : 
4? brs a? . a 


| ERO, Y. 


ca: un es 


= 
A = 
ry 
re Pe 
i 7 - 


| 5. oe RRS 
| , CRE = 0 tw | ‘ont 2 pf adtay tt omens 
rt a. te 


he 4 Linkin’ to wilerakebo 


se 


vo 


re ae 7 é if: a ; 
— emia BA. vs 08 Fuca at 2 


oe 


eu, 


solution. Therefore, besides the potential Vol Which is essen- 


tially repulsive near the origin, - V is also a possible po- 


Ol 
tential curve which is attractive near the origin. As an esti- 
mate of the quantum mechanical bound states with this attractive 


potential, - V was approximated by a square well potential: 


Ol 
Vo1 SP ie AG) 3G ae =. 

dimensionless units 
Voi = O se 2.5 


The minimum well depth for forming one bound state in a potential 
well of depth, V, and width, a, is given by: 
fp 
va2 >oe 


Converting to dimensionless units: 


2 
2 T G 
Vqaq > “3° 


Inserting numbers, the left hand side gives 8.74 while the right 
18° 16.7 ; V2.5 is too small by a factor of 2, which is signifi- 
cant even with the crude square well approximation, to give a 
quantum mechanical bound state. 

The total energy of the classically bound particle in 
this potential is found to be En = 5.47 me, with a binding ener- 
gy, B= 1.82 me®, 

Teshima also has calculated the potential energy curve 
for the interaction of two particles of the first kind, l.e., 

% and Po: This energy curve is called Voo and is reproduced 
from T, and shown in Figure 6. This potential is essentially 


repulsive near the origin, but if the sign of one of the solu- 


8S 


-eee0 et dottiw Vv v cetasieil aie wenn dd eetlanchall — 


~oq sftt#eoq’s offs el (AV = yfilgtto ort asen ovtetugen 
-Itas ns 8A .Afeiao elt TAR evitoatste af 


aVitosttis effi dttw aetetd prived Ixotnertoom my Sale a 
:igtinodoq [few ots¥pe s Yd bedamixomgus esw tov - pvt 


ive) a ae Ay % 
2.0 > %> 0 OF - = ro” ; 7 = ay 


y a 


@.0<¢ 1% 0 = so" 2 ; Dae 


ating aasinotenomrh | 


— 
- 


_Istineteq s at stata baved eno gntmrrot tot (qe om mutate At 
yd aevig et .6 .ddbiw bas Stqpb to Liew 


, ees aV ar 
“petkays avofacnartemtb eo : 
as we 


S 
ey 4 ; - a ad 
) < i ; res 


dogit ond oftnw SY,8 asvia abt. bid sitet oct yetidaum getasent 
-ttingte ef sotdw .S to todos s yd Lieme oot at fagh mar at 
s svin of .noty¢emixorqgs Iisw ovsupe sbvto aed _—— 


.ogate bimod: eotawdoon mud 
at elotéteq Dauod es one ‘to baie ee 
-rene antontd 6 déiw . om TH.8 = 48 9d 08: | 


_ 
= - —s 
See 
A 7 


a 
aed 

ae a 7 
el oe 


eee ; aad 


a ee pre c 
on _—. 


evino yEtene te at nat : 
pet ae fet Lotonaq owt 


Figure 6 


The interaction potential energy between two 


solutions of the first kind is given by the 


solid curve. Reproduced from Teshima, (1960). 


The dashed curve is the exponential approxi- 


mation. 
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tions %o is changed, the potential will be attractive. To in- 


vestigate the possibility of forming a quantum mechanical bound 


sravector ‘this situation, « Voo was approximated with the func- 
tion 
-r_/a 
as rr 
Ve=- Bs e ; 
Lay 


on -) a ee 
ag = 20) agree 
then V approximates - Voo quite well as is shown in Figure 6, 
The analysis of the Schrodinger equation for a potential of the 


ror, V= 5 ete Shows that no bound state will be formed un- 


less the condition 
CUB 
Re 


holds. The constant 2.405 is the first zero of the zero'th order 


eae} > 2.405 


bessel function of the first kind, and uw is the reduced mass. 
Using the dimensionless units together with the fact that the 
masses of the two systems - Po and % 5 are equal, the condition 
is: 5 
Baaq > Ee) " 

The numerical values for age Ba> and G give 13.3 for the left 
hand side and 19.5 for the right. Again the well is probably not 
deep enough to have a bound state. 

The binding energy for the classically bound particle 
is equal to the energy of the two separate parts and therefore 


the total energy is zero. 


The potential energy curves Voo and Vo1 still show the 
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possibility of forming several other classical bound states. 
These states occur in both cases for a separation of the particles 
of the order of 4X and have a very small binding energy. Since 
the potential curves used are only approximations, the existance 
of these bound states is doubtful and are listed only for com- 
pleteness. A more accurate investigation of the interaction po- 
tentials will have to be made to determine if these states exist. 
A summary of the results obtained above is shown in 
Tables 2 and 3. Further discussion of the physical interpretation 


of these bound states will be found in Chapter 6 of this work. 


Table 2 


Summary of Bound States 


Quantum Mechanical Classical 
Bound States Bound States 
per 
one 
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Table 3 


Doubtful Bound States 


Quantum Mechanical Claaeios) 
Bound States Bound States 
a 


none 


Interaction 


binding energy of system of particles composed of parti- 


to 
II 


ole 
cles 9,(- ®;) 


E = total energy of the system 
7? = classical separation of the particles 


B and HE are in units of fee with the coupling constant G = 13.53 
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IV. The Variational Method (V.M.) 


It is well known that the solution to a linear diffe- 
rential equation is the extremum of a certain integral, usually 
called the Lagrangian, where several subsidiary or auxiliary con- 
ditions may be required. In general the differential equation 
is the Euler equation for the Lagrangian. 

It is easily seen that equation (1-9) is the Euler equa- 


tion for the Lagrangian, 


2 2 1 4 
L = z{[- + (Vo) Oe ] dv’... (4-1) 
The integral is taken over all space and o must satisfy either 
o—O0O as r— ow, Or oe — 0 as r — >on. se is the normal de- 


rivative of » on a surface. Thus, the solutions of equation 
(1-9) satisfying the boundary condition (1-10) are the extremals 
of the above Lagrangian. 

It will be shown later that the constant, fo in (4=1) 
is chosen so that L will be equal to what has been defined as 
the energy of the field. 

The difficulty with the above variational method is 
that the second order variational terms are not positive or nega- 
tive definite, so one has no assurance that the method provides 
an upper or lower bound to L. This is shown in Appendix B. 

In order to remove the above difficulty concerning the 
bounds of L, it is necessary to reformulate the variational prin- 


ciple in a different way. Before doing this, it is necessary to 
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know that the solutions of equation (1-9) satisfy the following 


integral relationship: 


3| [v8 + (99)? ]av is Jotav 


This equation together with the following three integral rela- 


tionships are proved in Appendix A. 


3 | eav = | eSav (4-4) 
12 [ av Pe 


Ss) Je?av 


AS a consequence of the last two relationships, 


dv (4-5) 
J (v0)?av (4-6) 


note that the energy of the field 9, interacting with itself, as 


calculated from the definitions (3-1) and (3-2) is 


Ei = - {ov2pav ; 


Then by Green's theorem, 


B= i J] (v0)2av . (4-7) 
The same expression is obtained by using the integral relation- 
ships (4-5) and (4-6) in the Lagrangian (4-1). Thus for an ex- 
tremal, the Lagrangian is positive definite and is equal to the 
energy of the field. 
To reformulate the variational problem so that the 
eigensolutions always make the Lagrangian a minimum we must con- 


Sider either of the two problems: 
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be 


zs Extremize L =~ p-dv (4-8) 


subject to the auxiliary conditions, 


ds 12 [ av - J otav = 0 (4-9) 
and 
K = 3 | oav . f (ve) Fav = 8: (4-10) 
a 2 2 
poe Extremize L = +7— J + (Vo) Jo (4-11) 
subject to 


M = 3 | fF + (v9)*| av - [olay =O. (4-12) 


To show that the extremals of the above variation pro- 
blems satisfy equation (1-9) use the method of undetermined 


Lagrange multipliers. 


1 Extremize L' = L + AJ + uk (4-13) 
where A and w are the Lagrange multipliers to be determined. 
The Euler equation becomes 


uv°o + (3 + 12A + Bu) - QP =O . (4-14) 


Multiply (4-14) by 9, integrate, use Green's theorem and the 
auxiliary conditions (4-9) and (4-10) to get A = 1/4 and the 
Euler equation is: 


uvep + (6+ L)Q-ZOP=0. (4-15) 


In the same way that the integral relation (4-5) was shown to 


hold for solutions of equation (1-9), it can be shown that the 
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solutions of (4-15) satisfy 


ay = 8(6 + au) [ 9?av ; (4-16) 
But (4-16) is required simultaneously with (4-9), therefore 


I. = -3/2 and the Euler equation is identical to equation (1-9). 
The second problem is easier. 


het Extremize L' = L + 2M (4-17) 
where A is the Lagrange multiplier. The Euler equation is: 


(# + 3a)(~ - Vp) - ange =o. (4-18) 


Multiply (4-18) by mp, integrate, use Green's theorem and the 
auxiliary condition (4-12) to obtain } = 1/4. The Euler equa- 


tion (4-18) is identical to equation (1-9) with this value of 2. 


To see that the extremals which approach zero at infi- 
nity make the Lagrangian in the above problems a minimum, note 
that in both cases the Lagrangian (L) is a positive definite 
quantity. If the auxiliary conditions are not imposed, the se- 
cond order variational terms of L are positive definite. This 
can be seen by an analysis analogous to that given for the La- 
grangian (4-1) in Appendix B. Now, if any number of auxiliary 
conditions are imposed, they can only serve to raise the value 
of the minimum. See page 407, Courant and Hilbert, Methods of 


Mathematical Physics, Vol. l. 
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L an absolute minimum is 9 = O. 


To obtain the remaining solutions it is necessary to 
develop a variational method for "excited" states. To do this 
a connection is needed between the two solutions, ?5 and Pos of 
equation (1-9). This connection is furnished by the so called 


orthogonality relation , 


3 ae ce - 
fer 5 av = {93 9, dv, (4-19) 
which is proved in Appendix A. 


Now consider the variational problem: 


a oe Extremize L - dz | [23 + (99)? dv (4.20) 


3 [e+ veg)? |aw 
a [[s P5 - e2 ae 


where P+ is known and satisfies vo, = P45 ~ = Pi: 


subject to 


— 
Il 


(4-12) 


| 
a2, 
Mm 
Q. 
<j 
II 
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and 


(4-21) 


ra 
| 

II 
oO 


Whether or not we regard 7; = O as the first solution 
makes no material difference as the condition (4-21) is satis- 
fied identically and problem III then reduces to problem Il. 
Thus, the first non-zero solution of equation (1-9) is given by 
that extremal of problem II which now gives the lowest value for 


the Lagrangian not equal to zero. 


To obtain further solutions, consider variational pro- 
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blem III again, but with the requirement that ? be equal to the 
non-zero solution found above. This time the orthogonality re- 
lation will be important, since it is not identically satisfied. 


To determine the Euler equation consider: 


Bait, Extremize L' = L + AM + uN (4-22) 
where A and uw are the Lagrange multipliers. The function P41 is 
known and satisfies vo, =e - 92. The Euler equation is 

3 se ne eee ay 


To determine the constant uw, note that under the transformation 
P5— - Po; L and M remain invariant. Further, N is also inva- 
riant since it is equal to zero. Therefore the Euler equation 
satisfied by - Po must be identical to the equation satisfied by 


®>5- But from (4-23) with 5 replaced by - 9, there results, 


- (2 + 3A) (a5 - Vos) + 2r02 + (3950, - oF) = 0+ (4-24) 


Adding (4-24) and (4-23) we have 


2 2 
(305 zs 7) 1 =O. 


Therefore w = 0, since 1 + O and of te 305. 


To obtain the value for A proceed as in problem II'!. 


Again \ = 1/4 and the Euler equation is identical to equation 
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Thus formulated, it is seen that the above variational 
problems can be used to find approximate solutions of equation 
(1-9). It is only necessary to chose trial functions which sa- 
tisfy the required auxiliary conditions for the Lagrangian to 
assume a value at its minimum which is higher than the actual va- 
lue. In the remaining part of this work several examples will 
be presented, but it should be pointed out that the above varia- 
tional problems have not been applied systematically, since most 
of the variational theory was formulated after working out the 
results of many trial functions. 

To be more specific let us now describe the content of 
the variational method (V.M.). This method attempts to approxi- 
mate the solutions of equation (1-9) by substituting a trial 
function, Pp lay sAn+++-a43 FP), which must satisfy certain auxi- 
liary conditions, into a Lagrangian. The Lagrangian is then 
minimized or maximized with respect to the parameters of the trial 


function. Let these parameters be As» then 
Lala, ) = L( 9m) 


The condition that the integral be made an extremum leads to the 


Simultaneous equations, 


oly 
as =0; (1 =1,2,...-n) , (4-26) 
ae 
The solution of these n equations determine the optimum values 


of Os and therefore the best possible trial function of the form 


used. The particular form used initially for the trial function 
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is suggested by some a priori guesses as to the nature of the 


exact solution. 

The V.M. is a powerful tool with which we can approxi- 
mate the solutions of a non-linear differential equation which 
cannot be solved. One great advantage of the V.M. is that we 
now have a method to attack the solutions of equation (1-9) which 
are not spherically symmetric. One must, of course, show that 
the solutions of more complicated symmetry exist, but this pro- 


blem has generally been ignored in this work. 
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38. 
V. Calculations Using the V.M. 


In this section we shall present the results of the 
analysis using the V.M. for several types of solutions to equa- 
tion (1-9). Before presenting the results for the specific func- 


tions we shall discuss two points of general interest. 


ia Green's function method of improving the trial functions. 
A method of improving the trial functions can be devi- 


sed by writing equation (1-9) in the form 


Vo - @ = Anp(B) (1-9a) 
where 
o(#) = - + pF) . (1-9b) 


The solution of equation (1-9a) is given by 


og(#) = [aC )p(#)ae" ey 
where G(7,r!') is the Green's function for the operator (ve aD) 
and p(r) is assumed known. Now let 
si 1 3) 
p(f) = - tox Op (Ps ay----a), (5-2) 
then (Ts Qj ++++G,) should be a better approximation to the ei- 


gensolution of equation (1-9). It was shown in T. that Pa indeed 
is a better approximation for the nodless spherically symmetric 
solutions. Note, however, that we do not have any general proof 
that after n such operations implied by (5-1), Pq will approach 
the actual solution. In any case, if Pq contains some parameters, 


they can now be determined by the V.M. Another procedure is to 
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determine the parameters in Pn first by the V.M., insert the 
numerical values, then find Pq from equation (5-1). The diffe- 
rence between the energies obtained by these two methods is 


small if the trial functions are spherically symmetric. 


rae Auxiliary conditions and a specialized form of the Lagran- 
gian (4-1). 

Besides the extremization conditions on the Lagrangian 
of a trial function, one or more integral relationships as well 
as possible orthogonalization conditions must be satisfied by 
the trial function. Therefore, for any trial function, Pm there 
are three kinds of equations which may be required. 

n Extremization conditions 
oe = 0 (a eae yea) (5-3) 


Three possible integral relationships 


) 
12 | 9 av = on dv (5-4) 
3 (er av = [ (Ven) av (5-5) 
3 [on ani \ ep dv (5-6) 


Orthogonality conditions 
3 9,av = 3ay ‘Cas ee (5-7) 
Om V4 5 ae Pm bart — Ck <r er Tee 


Note that any of the auxiliary conditions used in the variation- 


al problems of the preceeding chapter may be formed from these 
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integral conditions. 

One of the extremization equations (5-3) is always equi- 
valent to the integral relation (4-12). This can be quickly veri- 
fied by letting Pn = An(®), where A is a constant, in the Lagran- 
gian (4-1). The equation ee = 0 yields the desired relation 
(4-12). Therefore the results of method II or III will apply 
if the above extremization has been used, even though the Lagran- 
gian (4-1) is minimized. 

A large number of the three parameter trial functions 
investigated reduce the Lagrangian (4-1) to a similar form. In 
this case, two of the extremization conditions (5-3) are equiva- 
lent to the two integral relations (5-4) and (5-5). To be more 
specific, the specialized form of the Lagrangian (4-1) will be 
obtained if the integrals in the Lagrangian reduce to: (The 


integrals (5-11) and (5-12) are included for future reference. ) 


Oo 
tz | % dv ae @, (p) (5-8) 
a ae a,(p) (5-9) 
ae |" oo 
be ne) od gees (p) (5-10) 
Hp | YOn/ OY = Bo MehP 
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te | % dv = 2 Q) (Pp) ane 


3 
3 3 Q, (P) (5-12) 


The Q,'s are functions of the parameter, p, but the restriction 
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to a single parameter is not necessary. B and 6 are the other 
two parameters. Note that B, 8, and p may not appear directly 

in the trial function (e.g. they may be combinations of the para- 
meters of the trial function, such as p = b/q.) The Lagrangian 


(4-1) will now be 


2 2 
L»(BsBsP) = 3 = au + 8°Q,(p) - > a(0)| . (5-13) 
The two minimization equations, 
) ) 
eaieas and ae ie (5-14) 
OB els) 
demand that the following two relations hold: 
a 
BT = 3Q,/@ (5-15) 
and 
2 
Bo = 12Q,/@ - (5-16) 


These values of 3° and a substituted into the special Lagran- 


gian (5-13) yield 
-— 2 ae (5-17) 


For a three parameter trial function, Ln is now a function of one 


parameter, the value of which is to be determined from the remain- 
ing extremization equation, 


re) 
= are (5-18) 


Alternately, a graph of the function Ln(p) can be made, the ex- 


tremums of L, giving the solutions of equation (5-18). Usually, 
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Ln is a complicated function of p, so that the graphical method 
is quicker and easier to apply. When the above method of solving 
for the parameters has been used, the designation (OR1) will be 
placed beside the energy thus obtained. 

If equations (5-15) and (5-16) are restored to their 
integral forms by using (5-8), (5-9), and (5-10), it is seen 
that the two integral relationships (5-4) and (5-5) are automa- 
tically satisfied. 

The only other method of solving for the parameters 
that is used here is called (OR2) and is described below. 

Suppose it is desired to have the integral condition 
(5-6) satisfied by the trial function. By using equation (5-11) 
and (5-12) it is seen that the integral condition may give a 
relation of the form, 


Bo = C,(p) . (5-19) 


A similar relation is always obtained from the orthogonality con- 
dition (5-7). This equation can be used in several ways together 
with the extremization conditions to determine the values of the 
parameters. The procedure (OR2) does the following. Substitute 


B° from (5-19) into the special Lagrangian, (5-13), reducing it to: 


C, (p) 
Ln(BsP) = 2 | B°@,(p) + Co(p) 
where (5-20) 
C, = @, - C1 &,/6 
3 


The extremization equation 3a- = O gives, 


g° =-3C,/@, - (5-21) 
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This value of ae used in (5-20) gives, 


In = 0,8, 305 (5-22) 


Ln is again a function of a single parameter which can be deter- 


mined as before. 


The results for various trial functions follow: 


dis Spherically symmetric solutions. 

The accuracy of the numerical solution Py Was discussed 
in T. In all of this work with the V.M. it is assumed that the 
numerical solutions Po and p, are more exact than any of the trial 
functions. Thus, there are three spherically symmetric solutions 
with which to compare the results of the V.M. calculations. In 
this work only two of these solutions have been approximated by 
the V.M. A table of the spherically symmetric trial functions 
appears below. The integrals and additional information for the 
functions appears in Appendix E. 

To show that the V.M. is giving reasonable numbers for 
the parameters, the functions 1S (OR1) and 9S (OR1) are plotted 
together with the numerical solutions Po and ?5 in Figures 7 and 
8, Notice that the V.M,. has determined the position of the node 
quite well for ?1 even when using the simple trial function 9S. 

The energy obtained by the V.M. is in all cases larger 
than the numerical value. This is to be expected for the solu- 


tions with no node where (OR1) has been applied, since (OR1) on 
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Table .4 


Spherically Symmetric Trial Functions 


Desig- Energy in 
nation Function units of mce“/G Remarks 


ie 


13.856 oo 
14.10 (OR2 


ea 
cylindrical 
coord tes 


ne ol 121 +p) 19.6 (OR1) ina 


18 
28 
Ts ” 
5S ae: 
68 | Aeo™ 4 ge: FE 14.07 (OR2) | on = V3; @f{0)= 0 
; if 
78 
ia : 
a 


13.526 “CORI see Fig. 
| = EP gives 
b= 0 


14.10 OR2 


?R 


ar 


89.745 Pon see Pig. 12 


A(1 - br)e™ 89.57 (OR2) | o, = WN3 


Note : The (OR ) following the energy refers to the method used 


for determining the parameters. See page 41 and 42. 
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The solid curve shows the first solution, Po» 
while the dashed curve is the approximation 


given by trial function 1S (ORL). 
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Figure 8 


The solid curve shows the second solution, 04; 
while the dashed curve is the approximate so- 


lution given by trial function 9S (OR1). 
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the Lagrangian (4-1) is equivalent to variational method I of 
Chapter IV. Note that (OR2) is not equivalent to any of the va- 
riational problems so that the energy given by this method may 
be above or below the exact value. The function 9S with one node 
also may not necessarily have definite bounds, since this func- 
tion has not been made orthogonal to the nodeless solution which 
is the necessary condition for the results of variational pro- 
blem III to apply. 

Function 8S is the one obtained by using function 1S 
as the source density in the Green's function method of improving 
the trial functions. Note that this function has the correct 
asymptotic behaviour for large r together with a vanishing deri- 
vative at the origin. Since the (OR1) energy for this function 
is the lowest attained by any of the trial functions in Table 4 
it is the "best" trial function. Table 5 shows a comparison of 
the numerical solution, 9), and the function 8S (OR1) for seve- 


ral values of the argument. 
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Table 5 


Values of the Function 8S (OR1) 
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Compared to the Numerical Solution Po 
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es Solutions of odd parity. 

We shall try to approximate the eigensolution of equa- 
tion (1-9) which is cylindrically symmetric, has odd parity, and 
only the one nodal surface required by its parity, assuming this 
solution to exist. Because of its parity, these solutions are 
automatically orthogonal to the spherically symmetric solutions 
so that orthogonalization is not required as a subsidiary condi- 
tion and variational method I or II may be applied. 

A table of the trial functions used appears below. 
Again, the integrals and parameter values for these functions are 
found in Appendix E. The table will not be discussed in detail 
Since it is self explanatory. 

Function 7P is obtained by using function 1P as the 
source density in the Green's function method of improving the 
trial functions. One would expect 7P to be quite a good approxi- 
mation to the actual solution, similar to the situation that oc- 
cured for the spherically symmetric solutions. However, the in- 
tegral relationship, (4-5), is not satisfied by a factor of about 
two. If the numerical work is correct, this probably indicates 
that the original source function, 1P, is a very poor approxima- 
tion to the actual solution. 

If the cos@ dependence of the trial functions is made 
more general, as in functions 6P and 8P, energies lower than the 
49,26 of function 1P are obtained. The only minimum for function 
8P is found for n—»«.. This seems to indicate that the field is 
concentrating itself in the x-y plane. 


Trial function 9P is the difference of two spherically 
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Tabile«§ 


Odd Parity Trial Functions 


Desig- 
nation Function 


Energy. in 
units of mc“/G Remarks 
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symmetric functions such as 1S separated by a distance, s. The 
Lagrangian as a function of 2as is shown in Figure 16. The func- 
tion L(s) has a minimum for s = », since qa approaches a constant 
for large s. Since the V.M. gives an upper bound, function 9P 
is the best solution of odd parity found so far. 

These results seem to indicate that equation (1-9) has 


no single particle solution of odd parity. 


cP Solutions with other types of spatial behaviour. 

Only a few functions with symmetries other than those 
previously discussed have been tried. These are listed in Table 
7 for reference only since most of them have not been investigated 


thoroughly. 
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Table 7 


Trial Functions with Various Symmetries 


Function 


a 
ee 


Desig- 
nation 


Energy in 
units of mc°/G 


Remarks 


Prolate spheroi- 
dal coordinates. 
c—>»0; reduces 

to 15S. 


13,01 OR) 


Oblate spheroi- 
dal coordinates. 
ec— 0; reduces 

CO: d3. 


1308. (ORL) 


c— 50; reduces 


' (@) f 
to y= Ar® (3r° 


118.24 (OR1) 
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Functions 1H and 2H are expressed in spheroidal coordi- 
nates and are therefore cylindrically symmetric. These functions 
were used in order to determine if the solution of lowest energy 
is spherically symmetric. The V.M. gives the reasonable result 
that the spheroid parameter, c, must be zero. Thus, functions 1H 
and 2H reduce to the nodeless spherically symmetric function 1s. 

Function 3H is another spheroidal type function which 
was used to try to find a solution with symmetry similar to the 
D state solutions of the hydrogen atom. The spheroid parameter 
again was found to go to zero with the spheroidal coordinates re- 


ducing to spherical coordinates. 
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VI. Conclusion and Discussion 


Since this non-linear theory is far from being compre- 
hensive at the present time, we should not expect the masses of 
the particles as calculated by the theory to compare very well 
with the experimental results. In fact they do not. However, 
we have been able to show how particle like solutions may be 
found which possess definite mass ratios of the same order of 
magnitude as those occuring for the elementary particles. By as- 
Suming an expression for the interaction energy, it is possible 
to form various bound states of these fundamental particles. The 
coupling constant, G, also assumes the reasonable value of 13.53 
if the energy of the first particle, Po» is taken to be me, 

The mathematical procedure of the variational method 
for the non-linear differential equation has been more rewarding. 
We have seen how the variational method may be used to find ap- 
proximate analytic solutions of a non-linear differential equa- 
tion. Furthermore, if the trial functions satisfy certain subsi- 
diary conditions, the minimized Lagrangian will always furnish 
an upper bound to the energy. By use of this variational method 
quite accurate analytic expressions have been obtained for the 
first two spherically symmetric eigensolutions. 

The variational method was used to try to find a solu- 
tion of equation (1-9) with odd parity. The trial function of 
odd parity, giving the lowest energy so far, is the difference 


between two nodeless spherically symmetric functions. The two 
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spherical functions must be separated by an infinite distance to 
give the lowest energy. In effect we have found that equation 


(1-9) seems to have no single particle solution of odd parity. 
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Appendix A: Proof of Integral Relationships 


In this appendix we shall prove the following integral 


relationships holding for solutions of equation (1-9): 


3foav = \ 3a (A1) 

3\ [9 + (¥e)* Jav = | ota : (A2) 
12 | 92av z f ate} (A3) 

3\ av = | (v9)%av , (a4) 


and if 9, and 9, are both solutions of equation (1-9), 

On = \ ogezav ; (A5) 
Integral relationship, (A3), was first noticed to hold for the 
numerical solutions of spherical symmetry and was proved analy- 
tically by Dr. H. Schiff for spherical symmetric solutions. It 
was thought it might fail for solutions of lower symmetry but 
Dr. D.D. Betts succeeded in proving it in general expressing the 
volume integrals in spherical coordinates. The following proof 


is carried out using the more compact vector notation. 


Proof ,of..(.Al): 
» satisfies equation (1-9) 


1 
V9 = 9 aie 9? (1-9) 
Integrate both sides of this equation over all space. By Gauss! 
theorem, Py 
V Qdv = " Vo ds 


V S 
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where the surface, S, bounds the volume, V. Since we consider 
only those solutions which are either zero or a constant at in- 
finity, the surface integral is zero if the volume integral is 


over all space. Therefore, 


which is (Al). 


Proof of (A2): 
Multiply (1-9) by » and integrate over all space. By 
using Green's theorem, 
ov & - | v0)2av ~ ic - 9! )av 
S 
If » is a constant or zero at infinity, the surface integral 


vanishes and equation (A2) follows. 


Proof of (A3) and (A4): 
Begin by integrating | 9av by parts; 
| otav =- 3 | (8 v)eav 
This relation can be proved by performing the integration by parts 
in rectangular coordinates. g must also be zero at infinity, for 


the non integral term from the integration by parts to vanish. 


Now use equation (1-9) to replace the o> term on the right giving, 
4 — ul geo(s 
mp dv = - 4 \ p(r-V)pdv + V“o(r-V) pdv 


After another integration by parts the first term on the right 


is 6 | aPav. The second term may be transformed to 
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= | Vo + V(¥-V)odv 
by Green's theorem. Using the vector identity, 
V(U-V) = U:VV¥ + VeVU + Ux(Vx¥) + Vx(Vxi) 
the integral becomes, 
- i vp. [ 00g + Vo-Vr + Px(VxVo) + Yex(Vx#) | av 


The last two terms are zero by virtue of the identities; 

VxVm = O and Vxr = 0 
The second integral is - 4 | (v9) 2av since the del operator on 
the radius vector, Yr, gives the unit dyadic. The first integral 
may be reduced by noting that 


Vo: (7%-VVo) = = PV(Vo 


2 
) 
This will give 6 | (ve) av after integrating by parts. Thus after 


collecting all the terms, 


| stay - 2 | (v9) av £ 6 | ePav (A6) 


If (A2) is now used the relations (A3) and (A4) result. Since 
Green's theorem as well as several integrations by parts have 
been used, p must be zero at infinity for the surface integrals 
and constant terms arising from the integrations by parts to be 


Zero. 


Proof of (A5): 
1 and 9, both satisfy equation (1-9) thus: 
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Vp = % - 3% (A8) 


Multiply (A7) by Pos (A8) by ®,, subtract and integrate over all 
space. Then by using the symmetric form of Green's theorem, the 
volume integrals containing the Laplacian reduce to surface in- 


tegrals which are zero if P+ and/or Po is zero or constant at 


infinity. The remaining integral is (A5). 
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Appendix B: Second Order Variational Terms 


The Lagrangian (4-1) for a trial function, > 1s, 
3 2 a 1 is 
ty = Be | [08 + (vom) merce ond 
Now suppose Pp = — + ef, where po is the actual solution of equa- 


tion (1-9), ¢ is a small parameter, and f is an arbitrary function. 


For this Pps the Lagrangian up to terms of second order becomes, 


Ln - 3; | [5% (Vp)* - é o'lav + 3 te = «fers + ev(T-Vo - = 5 ot |av 


at Es ee GES un 26°] av 


The first term on the right is L, the Lagrangian for the actual 
solution. By using the fact that » satisfies equation (1-9) and 


Green's theorem, the second term is zero. Thus Ln reduces to 


In - b= & | [es + : aye? | dv 


This integral is not positive or negative definite for arbitrary 
f as can be seen by chosing f = e7 and using the fact that 
{ g?av is bounded together with the Schwarz inequality. The sign 


of Ln - L will depend on the value of a. 
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Appendix C: Results of Improving Odd Parity Trial Function by 


Means of Green's Function 


let us find the solution to the equation 


Vo - 9 = -= 92 
where Qn = Arcosée %; A = eve ao a 
The solution to 

Vp - 9 = - Unp(r) 


ie given by 


p= | a(n.) p(eyav 


Fale = rt 
where G(r,r!) = ———-;y5\ k =i 


Jr - r?| 
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To perform the integration it is useful to use the expansion of 


G(r,r') in spherical coordinates: 


Where: 


P’ = Legendre functions of the first kind 
= spherical bessel functions 


e = 23 Mes 0 


Svge Mages 


For the definition of these functions see Morse and Feshbach, 
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Methods of Theoretical Physics. gis found in the form: 


2 
9 = -- 2, - yz) cose ~ p3cos7o| 


with 


II 


- 2x 10°y, = (5.7950889) (Z$+)e7_ (4 


diy a 6 


ye79 196r [76.9277 


S 122.99723r' am 117 .90271r° i 75 .33156r° + 30.111282r 


+ 5.7950889| 


3 2 
(0.335965) (F_+ SP + dor + 15) (=z) 
i 


6 s) 


a-oei96r 


: 
- 2x 10-y 
5 
i 


Mm 


| 76.927772" + 122,99720r° + 147492049? + 146.29556r 


+ 115.20744r> + 67.524951r° + 26.185185r + 5.039490 | 
A graph of the functions Vy and V3 ie Pound in Figure 9, 
To calculate the energy for the improved trial func- 

tions we have: 

vo a eels = 9 

Ji 

Multiply by 9, integrate over all space and use Green's theorem 
to give: 


\ccv9)? + p)av = + 3) 003 dv 


The energy is given by 


which can be reduced to 
i, 3 iT 4 
E - §;| 93 ay -t5\ ° oY 


by using the above relation. By numerical integration it was 


found: 
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Figure 9 


The functions y,; Yz, and (Wy - Y3) found by 
improving the odd parity function, 32rcosee 7% 3F 
by @ Green's function integration, . The curves 


Hix Wey and (Y4 - V3) are labled on the graph. 
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#;| eav = 49.452 


t=:\ p'av = 77.989 
By direct integration: 
te | 993 dv = 215.678 


Therefore the energy is: 


E = 107.839 - 77.989 = 29.85 
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Appendix D: Parameters for the Morse Potential Curves 


The parameters for the three Morse potential curves 


Shown in Figure 5 are given below. 
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Appendix E: Integrals of the Trial Functions 


With most of the trial functions used, the Lagrangian 
was reduced to a function of a single parameter with the func- 
tion being evaluated numerically to determine the extrema. There- 
fore, most of the results are Significant only to three figure 
accuracy. If the extrema have been found analytically the fact 
is stated in the "Results". 

To simplify the notation the following abbreviations 


have been used: 


P = tz \ ofav 
Q = | oav 


af 
a= tz Pndv 
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Se ire | on8y 


i 
L=P+R - ¢Q= L,(o,) 


E = value of L at minimum or maximum 
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Results: 
(OR1) 
E = 24AN3 = 13.856 
Am = 96 
a = 3 
+- 
E = 61AN33 = 14,10 
A = 9.0 
a = N33/4 
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